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Abstract: We consider a model for neutralino dark matter candidates arising 
in anomaly-mediated supersymmetry breaking schemes, and examine its testability 
through the search for exotic cosmic rays produced by neutralino pair annihilations in 
the dark halo of the Galaxy. We find that the model is already constrained by avail- 
able antiprotons and positrons measurements and may be further tested in upcoming 
measurements of these cosmic ray species. We show also that the monochromatic 
gamma-ray flux from neutralino annihilations is enhanced in this model up to two 
orders of magnitude with respect to alternative scenarios. The gamma-ray flux de- 
tected by the Energetic Gamma Ray Experiment Telescope in the direction of the 
Galactic center exceeds significantly the theoretical expectation of standard emission 
models. We prove that if at least 10% of this excess is due to gamma-ray radiation 
with continuum energy spectrum from neutralino annihilations in the model under in- 
vestigation, the associated gamma-ray line will be detected by upcoming gamma-ray 
experiments. 
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1. Introduction 

The low energy phenomenology of a supersymmetric version of the standard model 
depends critically on the mechanism to describe supersymmetry (SUSY) breaking, 
still an open question. In a recently proposed scenario, SUSY is broken dynamically 
in a hidden sector with dominant SUSY breaking terms induced by anomalies [|], 
|2j. Generic predictions of the anomaly mediated SUSY breaking (AMSB) scheme 
are the presence of a very heavy gravitino, with mass naturally in the 100 TeV 
range or heavier, and the fact that the lightest neutralino, plausibly the lightest 
supersymmetric particle (LSP), can be either a Wino or a Higgsino. 

One of the byproducts of SUSY is the possibility for the LSP to be a dark 
matter candidate. From this point of view, the case of Wino- or Higgsino-like LSPs is 
generally regarded as not very attractive: Their thermal relic abundance is typically 
rather small; it is in the range compatible with the observed dark matter density 
only for very massive neutralinos, as heavy as 1 TeV for Higgsinos, see, e.g., Ref [^], 
and a few TeVs for Winos, see, e.g., Ref. M. In the AMSB scenario, however, the 
nowadays density of the LSP is not fixed by its thermal relic density: An additional 
"non-thermal" source is present due to decays into LSPs of gravitinos or moduli 
fields, fields that parameterize a flat direction of the theory and that dominate the 
energy density in the early Universe [§, [J. Fine-tuning in this mechanism is needed 
for the LSP to be the main dark matter component, certainly a drawback of the 
AMSB scheme. On the other hand, we wish to stress here that the scenario with 
these LSPs accounting for the dark halo of our Galaxy has several phenomenological 
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implications and will be probed by upcoming experiments for large regions in the 
SUSY parameter space. 1 

There are several complementary techniques to identify dark matter neutralinos 
(for comprehensive reviews, see Refs. |7|, §). Direct detection of AMSB candidates 
was discussed in Refs. || 10 1. We focus here on another technique: the search 



for exotic cosmic rays, such as high-energy gamma-rays, antiprotons and positrons, 
produced by neutralino pair annihilations in the Galactic halo [IT|. This method 
looks particularly promising for neutralino dark matter candidates in the AMSB 
scenario, as both Winos and Higgsinos have a very high annihilation rate into gauge 
bosons, giving rise to strong sources of exotic cosmic rays. Indirect detection of 
Wino dark matter through antiproton and positron cosmic ray measurements was 
mentioned in Ref. ||; here, in a generic AMSB scenario, we derive the constraints 
on the model from available data and comment on the detection prospects for future 
experiments. Even more promising is indirect detection with the next generation of 
7-ray telescopes. We compute the monochromatic gamma-ray flux [12|] from dark 
matter neutralino annihilations, a signal with no conceivable background from known 
astrophysical sources, and show that it is greatly enhanced for AMSB models. We 
compare with present data and en-light the cases in which such flux will be detected. 

The outline of the paper is as follows: In Section ^| we describe briefly the particle 
physics model we consider. In Section |3] we discuss implications of antiprotons and 
positron fluxes, while in Section f| we consider gamma-ray fluxes. Section |] concludes. 



2. The MSSM in the AMSB scenario 

We work in the minimal supersymmetric standard model (MSSM) as defined in 
Refs. |~3|, |7j; details on our notation are given in Ref. ||, while, for quantitative 



prediction, we use the DarkSUSY computer code |14j. We suppose the LSP is the 
lightest neutralino, defined as: 

Xi = N n B + N 12 W 3 + N 13 H° + N 1A H° 2 . (2.1) 

The coefficients Nij, and hence the gaugino fraction Z g = |-/V n | 2 + |A 12 | 2 , are obtained 
by diagonalizing the neutralino mass matrix; they are mainly a function of the bino 
and the wino mass parameters Mi and M 2 , and of the Higgsino parameter [i. If \p\ 
is much smaller than \M\\ and |M 2 |, the lightest neutralino is Higgsino-like and Z g 
is close to zero. In the AMSB scenario the gaugino mass parameters are predicted 
in terms of the gravitino mass and their ratio is set by the corresponding gauge 
coupling constants and /3-functions coefficients of the gauge coupling constants. At 
the electroweak scale this relation becomes Mi ~ 3 M 2 , an assumption we keep 



lr The model has several more virtues; e.g., it solves the so-called "cosmological moduli problem" , 
see Ref. 0. 
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throughout the paper. If |M 2 | is much smaller than \p\, Z g is close to 1 and the 
lightest neutralino is Wino-like. We recall that, on the contrary, in alternative SUSY 
breaking scenarios, the GUT relation M\ ~ 0.5 M 2 is usually assumed and a gaugino- 
like lightest neutralino is forced to be Bino-like. 

For our purpose, it is not necessary to specify the whole mass spectrum according 
to a specify AMSB model. The cosmic ray yield per LPS pair annihilation is fixed 
by the annihilation cross section and its branching ratios. For both Winos and 
Higgsinos, the annihilation rate is fairly large (total rate av ~ 10 _24 cm 3 s _1 ), and 
dominated by gauge boson final states. Our predictions are then insensitive to most 
parameters in the sfermion and Higgs sectors, unless the annihilation cross section is 
enhanced further in some special configuration. We make a sample parameter choice 
which removes this possibility: We assume that all sfermions are degenerate in mass, 
with masses 10 times the lightest neutralino mass and with no-mixing between right- 
and left-handed components. In the Higgs sector, we fix the pseudoscalar Higgs A 
to be heavy, uia = 500 GeV and keep tan/3 as a free parameter. Finally, to specify 
a model, rather than assigning M 2 and /i, we fix the lightest neutralino mass and 
gaugino fraction M x and Z g , and the sign of \i. 



3. Antiproton and positron fluxes 

Although antimatter seems to be scarce in the observed Universe, a small flux of 
cosmic ray antiprotons and positrons is expected from the interaction of primary 
cosmic rays with the interstellar medium (see, e.g., Ref. [15]). The measured p flux 



[T7| and references therein) is compatible with the standard prediction for this 
secondary component, while strong exotic p sources can be ruled out. An analogous 
conclusion follows from e + data ( ||18|| and references therein). 

The neutralino-induced source of exotic cosmic rays is not fully specified once 
the particle physics model is chosen. The source function is proportional to the 
square of the neutralino density locally in space, a distribution which is only loosely 
known. In case of charged cosmic rays, the induced flux at earth is dominated by 
nearby sources; hence, for a given AMSB model, constraints on the local distribution 
of neutralinos can be derived. 

We estimate the neutralino induced p flux according to the analysis in Ref. []HJ, 
in which an accurate simulation of the antiproton yield per annihilation is performed 
and the propagation of antiprotons in the Galaxy is treated in a diffusive two- zone 
model. We suppose first that neutralinos are smoothly distributed in a spherical 
dark halo and compare with the data collected by the Bess [0 and Caprice |L7| 



experiments in their latest flight to derive the upper bound p on the local neutralino 
density. In Fig. |I] we plot, with a solid line, a few isolevel curves for p in the 
plane gaugino over higgsino fraction Z g /(1 — Z g ) versus neutralino mass M x (fi > 
and tan/? = 3 are assumed; alternative choices give hardly distinguishable curves). 
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Figure 1: Isolevel curves for the maximal local halo density p (solid lines) and maximal 
dumpiness parameter f5 (dashed lines), for neutralino dark matter candidates in the 
AMSB scenario, and in the plane gaugino over higgsino fraction versus neutralino mass. 
Limits for lower neutralino masses (on the left-hand side with respect to the dotted line) 
are derived from antiproton cosmic ray data, while those for larger neutralino masses follow 
from positron cosmic ray data. 



The displayed limits are at 90% CL. and are derived with a x 2 method under a 
few conservative assumptions. We compute the neutralino-induced p flux in all the 
energy bins where data are available, but include in the statistical analysis only 
those bins in which the expected flux exceeds the measured flux; in such bins the 
secondary contribution is supposed to be subdominant and is neglected. The values 
of p displayed are derived assuming the dark halo is described by an isothermal sphere 
with a 5 kpc core radius; a density profile singular towards the Galactic center would 
give more stringent limits []19~| . Finally, to model the diffusion of antiprotons in 
the Galaxy, we choose a standard set of the parameters in the propagation model 
compatible with measurements of other cosmic ray species (default choice in ||f 9|| ); 
extreme values for these parameters may enhance (as well as suppress) the flux up 
to about 50%. If one of these three assumption is changed, inducing a scaling of 
the maximal allowed p flux by roughly a factor F, the values assigned to the isolevel 
curves should be scaled by \f~F~. 

Fig. [I] shows that, for a given neutralino mass, Wino-like neutralinos (at the top 
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in the figure) have more stringent limits than Higgsino-like (at the bottom) or mixed 
neutralinos. The lowest upper bound on the local density is about 0.06 GeV cm -3 
for a pure Wino with mass about 100 GeV. There is a large region, with M x between 
80 and 200 GeV, in which p is lower than the value for the local dark matter density 
inferred from dynamical measurements in the Galaxy, about 0.3-0.5 GeV cm -3 if 
the dark halo is assumed to be spherical P0" |; an additional dark matter component 
is then needed. 



N-body simulations of hierarchical clustering |2T| suggest that dark matter may 
form substructures, or "clumps" , rather than being smoothly distributed in the halo 
(the picture we implemented so far). Larger local neutralino densities would enhance 
the sources of exotic cosmic rays [p^] . Suppose that a fraction / of the dark matter 



forms clumps with a some typical overdensity 5, defined as: 

s = 1 Jd 3 r d (pdjrd)) 2 

0.3 GeV cm" 3 / d 3 r d p c/ (r^) 



(3.1) 



where the integration is over the extent of the clump and p c \ is the neutralino density 
profile in the clump (see [f^| for further details). In this case, the induced cosmic-ray 



fluxes scale with the product / 5: comparing with the data, we have derived maximal 
values for this quantity. Isolevel curves for the maximal dumpiness parameter fS 
are shown in Fig. [I] as dashed lines. For heavier and heavier neutralinos the limit 
inferred from the p flux becomes less and less stringent. We find that, for M x between 
400 and 600 GeV, on the right hand side with respect to the dotted curve plotted in 
Z g /(1 — Z g )-M x plane, the analogous limit from e + cosmic ray measurements is more 
restrictive. This limit is derived computing the e + flux according to the analysis in 
Ref. |^| and comparing it with data from the Heat experiment [18], with the same 



approach as for the p flux. 

In the next years, space-based experiments |25, are going to measure the p 
and e + fluxes with better statistics and much wider energy coverage. It will then be 
possible to set even more stringent limits; furthermore, data at high energies may 
give evidence for the presence of this exotic component. As mentioned, AMSB mod- 
els annihilate mainly into gauge bosons and these channels automatically produce 
distinctive features in the high energy spectra if the neutralino-induced fluxes are at 
the level of the secondary components: a bump in the e + spectrum is expected at an 



energy roughly equal to M x /2 |27], |24|, a break in the p spectrum may appear at few 



tens of GeV |28[ . In addition, a low energy exotic component in the p flux might be 



identified [O. 



4. Gamma-ray flux 

Gamma-rays with continuum energy spectrum are generated by the decay of 7r° 
mesons produced in jets from neutralino annihilations. It is unfortunately a rather 
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featureless flux, difficult to discriminate from other plausible sources. When com- 
pared to the measured high-latitude gamma-ray background, it gives limits on the 
neutralino distribution in the AMSB scenario which are less restrictive than the 
correspondent bounds we derived from p and e + measurements. A much better sig- 
nature than the continuum contribution is given by the monochromatic 7-ray lines 
which arise from the 1-loop-induced neutralino annihilation into the 27 and Z'j final 
states (the energy of the 7-rays in the final state is equal, respectively, to M x and 
M x — Mz/4M X ). Only recently, the amplitude of these processes was computed at 
full one loop level and for a generic MSSM [|9|, 0. Predictions for these two cross 
sections, in the AMSB scenario and as a function of neutralino mass, are shown in 
Fig. H]. We have chosen a few values of Z g (one for each shaded indicated 
in the figure) and varied tan/3 between 1 and 60 (spread in each shaded area). The 
main contribution to the cross section comes from diagrams with charginos and W 
bosons in the internal loops; these are sensitive to the mass splitting between the 
lightest neutralino and the lightest chargino, which is slightly modified by tan (3. In 
all cases the cross sections are fairly large; for Wino-like neutralinos they are a cou- 
ple of orders of magnitude larger than the maximal values in other SUSY breaking 
scenarios |H|]. The case Z g = 10~ 2 corresponds to a good accuracy to the lower limit 



for each of the two cross sections; this is again a peculiarity of the AMSB scenario, 
as, on the other hand, there is essentially no lower bound to both cross sections for 
Bino-like dark matter candidates. Note also that in most cases the annihilation rate 
into the Z7 final state is larger than twice the annihilation rate into two photons. 

The 7-ray flux in a given direction of the sky is obtained summing the contri- 
butions along the line of sight; the largest flux is expected in the direction of the 
Galactic center, where an enhancement in the dark matter density is usually pos- 
tulated. The distribution of dark matter in the inner part of the Galaxy is still a 
controversial issue. N-body simulations find singular profiles, scaling as 1/r [32] or 
1/r 15 [33] as the galactocentric distance r — > 0. These profiles correspond to snap- 



shots of the Galaxy before the baryon infall; the appearance of the massive black 
hole at the Galactic center and of the stellar components may have sensibly modified 
such pictures with further enhancements (but a depletion is possible as well) in the 
central dark matter density ^4] |3lJ . There may be even the possibility that density 
of neutralinos in the galactic center region is substantial, while their contribution to 
the local dark matter density is subdominant. We parameterize the dependence of 
the 7-ray flux on the neutralino distribution through the dimensionless quantity: 

8.5 kpc Ji.o.s. (0.3GeVcm 6 y 

where p(l) is the neutralino density at the distance 1 from the Earth along the line 
of sight. Then , we introduce j = (J (tp — 0))&n = i -3 SI , the average of J over the 
field of view AQ = 10~ 3 sr and in the Galactic center direction ■0 = 0. E.g., the 
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Neutralino Mass [GeV] 



Figure 2: Annihilation rate into the two photon and the photon-Z boson final states for 
neutralinos in the AMSB scenario versus neutralino mass; a few values of the gaugino frac- 
tion Z g are considered. Also shown are 5a sensitivity curves for a generic next generation 
air Cherenkov telescope p6| , |37| , |38| observing the Galactic center: Dotted lines correspond 
to a few sample dark matter distribution, including a Moore et al. profile |33|] for j = 10 5 
and a NFW profile 32] for j = 10 3 . Dashed lines are obtained, for each value of Z g , by 
fixing the neutralino distribution in such way that the neutralino-induced 7-ray flux with 
continuum energy spectrum is compatible with the excess found in Egret data [41|. 



Moore et al. profile scaling as 1/r 1 - 5 and normalized to 0.3 GeV cm -3 at the 
Sun's galactocentric distance, gives j equal to about 10 5 ; for the Navarro Frenk and 
White profile scaling as 1/r and with the same normalization, j ~ 10 3 . For a 
given j, we can derive the sensitivity curve for a typical air Cherenkov gamma-ray 



telescope presently under development 



we consider an instrument with 



a 10 9 cm 2 effective area and 15% energy resolution, assume a 10 6 s exposure, and use 
a standard estimate for the background (see ]31j for details). We plot the relative ha 
sensitivity curves for a few values of j in Fig. § as dotted lines; sensitivity curves for 
the upcoming Gamma-ray Large Area Space Telescope (Glast) are comparable. 
As it can be seen, even for moderate values of j, large portion of the parameter space 
will be tested. 

The Energetic Gamma Ray Experiment Telescope (Egret) has mapped the 
diffuse Galactic gamma-ray flux up to an energy of about 20 GeV EI. The flux 
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detected in the direction of the Galactic center exceeds significantly the theoretical 
expectation of standard 7-ray emission models One of the conceivable interpre- 
tation of the data is the presence of an exotic component from the annihilation of 
dark matter particles |^T[. If we invoke the case of AMSB models, for given M x and 
Z g , we can infer the maximal value of j such that the neutralino- induced gamma-ray 
flux with continuum energy spectrum is consistent with or do not exceed the flux 
inferred from the data analysis. The line and continuum energy fluxes scale accord- 
ing to the same factor j; we can now consider this upper limit, say j max (M x , Z g ), 
derived from the continuum energy flux, and compute the corresponding sensitivity 
curves for the monochromatic fluxes (same assumptions as in the example above). 
These curves are shown in Fig. |2| as dashed lines, one for each value of Z g . Com- 
paring each curve with the prediction for the line cross sections, we can formulate a 
rather strong statement: If at least 10% of the excess found by Egret is due to the 
gamma-ray radiation with continuum energy spectrum from a neutralino dark mat- 
ter candidate in a AMSB model, the associated gamma-ray lines will be detected by 
future gamma-ray experiments for any values of the SUSY parameters in the model. 
Finally, by comparing the sensitivity curves obtained for a fixed j with those ob- 
tained by fixing the continuum energy flux, we find that some combination of SUSY 
parameters and dark matter distributions are already inconsistent with gamma-ray 
measurements in the Galactic center direction. Consider one sensitivity curve for a 
given j (dotted line) and one labelled by a value of Z g (dashed line): values of M x 
between the intersection points of the two curves are excluded (unless we invoke an 
efficient mechanism to absorb GeV gamma-rays in the Galactic center region). 



5. Conclusion 

We considered models for neutralino dark matter candidates arising in anomaly me- 
diated supersymmetry breaking scenarios, and discussed their testability through the 
search of exotic cosmic rays produced by neutralino annihilations in Galactic halo. 
We have shown that available cosmic ray data already place significant constraints, 
while even tighter limits are expected with data from upcoming experiments. At 
the same time, new data, both on the antiproton and positron cosmic ray flux as 
well as from gamma-ray surveys, may lead to the identification of these dark matter 
candidates. We proposed indirect detection through the search of a monochromatic 
7-ray flux in the direction of the Galactic center. We showed that such flux may 
be observable if there is a moderate enhancement of the dark matter density in the 
Galactic center region, and will be detected if the excess in the gamma-ray radiation 
from the Galactic center found in Egret data is caused by self annihilations of dark 
matter candidates in this class of models. 
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